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Abstract. Energy dependence of quark transverse flow carries information 
about dynamical properties (equation of state, initial conditions) of deconfined 
matter produced in heavy ion collisions. We assume quark- antiquark matter 
formation in Pb+Pb collisions at CERN SPS and Au+Au collisions at RHIC 
energies and determine quark transverse flow at the critical temperature of the 
quark-hadron phase transition. Coalescence of massive quarks is calculated in 
the MICOR hadronization model and hadronic final state effects are considered 
using the GROMIT cascade program. Comparing theoretical results to data, 
transverse flow values are determined and energy dependence is discussed. 
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1. Introduction 

In heavy ion collisions the formation of quark-gluon plasma (QGP) is expected 
at RHIC energies and above, however at lower energies we may not see a fully 
deconfined state. Although the precursors of the QGP state has been seen at 
CERN SPS [llj, it is still controversial whether this new state of matter could really 
appear in this energy range. To retrieve information about the energy dependence of 
quark matter formation, experimental data must be compared to model calculations 
based on deconfinement in a wide energy range. Experimental data on transverse 
momentum spectra for different particle species existed at ^snn — 17.3 GeV [E], 
130 GeV [El m and 200 GeV [0 0. Recently, NA49 Collaboration published new 
data at energies ^siS" = 8.8 GeV and 12.3 GeV {Eu^ ^AOA GeV and 80 A GeV) [ 
ISl El ) which extended the accessible region. 



1219-7580/04/ $ 20.00 
© 2004 Akademiai Kiado, Budapest 



2 



P. Csizmadia and P. Levai 



In this paper we are interested in the transverse dynamics of heavy ion coUi- 
sions. Transverse collective flow is directly related to the pressure of the central 
region with high energy density. Thus any jump (or discontinuity) in the flow can 
indicate modifications in the equation of state. We investigate hadronic transverse 
momentum spectra in the energy range ^snn — 8.8 — 200 GeV and determine quark 
transverse flow before hadronization. 

We assume the formation of deconfined matter with high energy density at 
all bombarding energies in the above mentioned energy range. Introducing critical 
temperature Tc for the quark-hadron phase transition, we assume a thermally equili- 
brated quark matter. The analysis of lattice-QCD results indicates the formation of 
a (constituent) quark-antiquark dominated deconfined matter at phase transition [ 
. This quark-antiquark matter can be described by macroscopical parameters as 
temperature, particle densities and flow velocities. Here we consider to be a uni- 
form phase transition temperature at all collision energies. Quark densities and flow 
velocities are energy dependent. Since quark flow cannot be measured directly, we 
determine it indirectly from hadronic momentum spectra. We apply hadronization 
models based on coalescence of massive quarks and antiquarks with small relative 
momentum (ALCOR model [ El E] to reproduce particle yields and MICOR [ 
1141 115j to reproduce transverse momentum spectra at low momenta, pt < 2 GeV). 
Our analysis is based on the MICOR model implementing different quark flow ve- 
locities and calculated results are compared to measured data. In parallel we will 
investigate the influence of hadronic rescattering on hadronic flow by the GROMIT 
cascade program [ITC]. 

Recently parton coalescence models were applied successfully at RHIC energies 
in the intermediate momentum region (2 < < 6 GeV) to explain the measured 
anomalous p'^ /tt^ and p~ /'!t~ ratios and reproduce the measured transverse mo- 
mentum spectra [El El El- In parallel, application of parton coalescence can 
reproduce the measured properties of the elliptic flow at RHIC energies [El El 1^ • 
The success of these calculations supports the application of quark coalescence in 
a wide momentum region. The MICOR model offers the possibility to investigate 
coalescence, especially flow phenomena, in low-px region. 

2. Calculations with the MICOR model 

In our model the expanding quark-antiquark matter cools down to critical temper- 
ature, Tc = 170 MeV, where hadronization starts. The longitudinally expanding 
deconfined matter is limited to the coordinate rapidity region rj = [— r^max, +?7max]- 
Inside this region we assume constant (anti)quark densities, dN^j/drj. We use the 
following notation for the numbers of newly produced quark-antiquark pairs at Tc- 

In the MICOR model final hadron numbers are scaled linearly on the initial 
quark numbers. This linear scaling remains valid for the rapidity densities also. 
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Fig. 1. Transverse momentum spectra for tt^ , , and p at y^SNN = 200 GeV 
in AuAu collisions (dots, squares and triangulars) [Ej- Full lines indicate quark 
coalescence results at vt = 0.60, neglecting final state interactions. Thick dashed 
lines show the applied exponential fits to extract slope parameters. Thin dotted lines 
indicate the MICOR results at a smaller transverse flow, vt = 0.55, but including 
flnal state interactions (FSI). 

In the ALCOR calculations for mid-rapidity data one can see the collision energy 
dependence of the number of newly produced quark-antiquark numbers [ I13| . which 
can be described approximately as dNqq/dy = 130 • ln(y^i]v7v/2.5) [ 1^ . 

The MICOR model covers wide space-time rapidity region and it reproduces 
measured mid-rapidity momentum spectra by the following parameters: ?7max — 1-9, 
Nqq = 720, fs = 0.35 at ./IWn = 17.3 GeV (SPS), and jy^ax = 4, Nqq = 3800, 
fs = 0.26 at ./JWn = 200 GeV (RHIC). At different SPS energies the above Nqq 
value is scaled by the measured charged particle yield, as well as at RHIC energies. 
(A detailed analysis will be published elsewhere [I24|.') 

The inverse slopes for transverse hadronic momentum spectra depend on the 
quark flow but not on the quark number. Secondary interactions among hadrons 
depend on absolute hadron numbers. We can set all above mentioned parameters 
of MICOR, but left only one parameter free, namely the quark transverse flow, 
Vt- We investigate the vt dependence of hadronic transverse momentum spectra 
of different hadronic species. In this comparison, as one of the possibilities, we will 
assume fast hadronization and neglect final state interactions. 

We consider 7r+ , , and p transverse momentum spectra at ^/snn = 200 GeV 
[El in AuAu collisions, see Figure ^ We expect an energetic collision with a fast 
expansion and fast hadronization. In this case we can compare the MICOR results 
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Fig. 2. Inverse slopes from MICOR model and RHIC experiment. MICOR points 
are connected by thin lines for vt = 0.55 ± 0.02 flow velocities and thick lines for 
VT = 0.60 ± 0.02. ./IWn = 130 GeV data points from PHENIX [E) are denoted by 
open squares, 200 GeV points from PHENIX [[^ are denoted by dots. 

(full lines) to the experimental data directly. The agreement between MICOR 
{vt = 0.6) and the data points at 0.2 < px < 1.5 GeV indicates the validity of the 
above expectation. Dashed lines show exponential fits to the calculated spectra. 

On the other hand, we are able to switch on hadronic interactions in the final 
state by using the GROMIT transport code ["TB*. However, GROMIT needs space- 
time information, which was integrated out in the MICOR model. Here a cylindrical 
symmetric matter was assumed with uniform densities, longitudinal Bjorken flow 
and constant transverse flow. 

The number of secondary collisions suffered by a particle moving with velocity 
V can be determined using the following rough estimate: 

oo oo 

n{t)(jvdt - I — avdt = , (2) 

J t"Wflow Wflow 
to to 

where n(t) is the density of the matter at the location of the particle at time t and 
a denotes the averaged cross section. The density, n{t), is supposed to be inversely 
proportional to the average flow velocity Uflow and a power of the elapsed time i". If 
the exponent a is a constant, then C is also constant. The average flow velocity in 
our case can be approximated by the transverse quark flow velocity vt- The result 
of this simple calculation is that a particle suffers at least one collision {Ncou{v) > 1) 
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Fig. 3. Transverse momentum spectra for tt"*", , and p at ^/snn = 17.3 GeV 
in PbPb collisions (dots, squares and triangulars) [13 |H|- Full lines indicate quark 
coalescence results at vt — 0.60, neglecting final state interactions. 



if its individual velocity v satisfies the following relation: 

V > /Co-. (3) 

Now switching on final state interactions (thin dotted lines) and focusing on 
proton spectra, we are able to reproduce proton data with a smaller quark flow 
superposed with the gained transverse energy during secondary interactions. On 
the other hand, the calculated light meson spectra will deviate from data, because 
the "pion wind" decreases the pion yield at higher momenta. From Figure ^ we 
can claim that there is no free room for final state interactions. Assuming fast 
hadronization, we are able to reproduce experimental data. However, if we want to 
include light mesons produced via independent jet-fragmentation, then the second 
case with final state interaction looks more favourable, where extra meson yields can 
be accommodated. (Contribution from jet fragmentation appears at much higher 
Pt for protons.) At this point we can not decide between the two investigated 
scenarios. Further analysis of perturbative QCD contributions is needed, similarly 
to Refs. [ini CHI Cn] and Ref. [E2- 
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Fig. 4. Inverse slopes from MICOR model without final state interactions and 
experimental data from NA49 [ [3 El QI| ■ MICOR points are connected by 
dotted lines for vt = 0.42 flow velocity, thin lines for vt = 0.44 and thick lines for 
vt = 0.49. Experimental data points at 40^ GeV are denoted by open triangles, 
80 A GeV data are displayed by squares, 158 A GeV data are indicated by filled dots. 

Figure 121 summarizes our results on slope parameters obtained in MICOR, when 
final state interactions are neglected. The uncertainties in quark transverse flow 
are displayed selecting vt = 0.55 ± 0.02 at ^s^n = 130 GeV (thin lines) and 
VT = 0.60 ± 0.02 at y^s^N = 200 GeV (thick hues). The MICOR results agree 
with the measured slopes quite well, which agreement confirms the application of 
MICOR model and the assumption of quark-antiquark matter formation at RHIC. 

At SPS we expect slower expansion and hadronization, which yield to a less 
energetic collision dynamics. Thus after hadronization, secondary interactions may 
have a larger role to reshape the final momentum distributions and modify the 
inverse slopes. To validate this expectation, we repeat our MICOR calculation at 
SPS energies with and without final state interactions. 

Figure (HI displays the MICOR result for 7r+, K^, and p transverse momentum 
spectra at ^/s^n = 17.3 GeV [[7\ |S] in PbPb collisions at SPS. Without final state 
interactions (full lines) one can reproduce the pion and proton data at vt = 0.49. 
For kaons we find slight deviation. Introducing final state interactions and focusing 
on the proton spectra, we obtain vt = 0.45. Proton spectra is reproduced in both 
cases. Kaon and pion spectra are similar where data are available. However, final 
state interactions modify pion and kaon spectra by decreasing the yield at higher 
Pt- The obtained suppression in the spectra at SPS is much less than at RHIC. 



Energy dependence of transverse flow in heavy ion collisions 



7 



400 r ' 1 ' 1 ' 1 ■ 1 ' 1 ' 1 ' 1 ' 1 

MICOR+FSI Vt=0.45 Vsnn=8.8 GeV, 1 2.3 GeV, 1 7.3 GeV 

375 - NA49, 158AGeV i » i 
MICOR+FSI Vt=0.41 



NA49, 80A GeV i— & 
MICOR+FSI Vj=0.38 




150 - JI K* p * A E" n' - 
'. 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 : 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 

m [GeV] 

Fig. 5. Inverse slopes from the MICOR model including final state interactions on 
NA49 experimental data [Ej |H1 (HI 1^1]' MICOR points are connected by dotted lines 
for VT = 0.38 flow velocity, thin lines for vt — 0.40 and thick lines for vt — 0.44. 
40 A GeV experimental data points are denoted by open triangles, 80 A GeV points 
are open squares, 158 A GeV points are filled circles. 



The displayed differences between data and the MICOR calculation with final 
state interaction indicate that hadronic collisions can become important in a wide 
energy range. However, we need to investigate whether pQCD processes (namely 
independent jet fragmentation) could provide the missing contributions both at 
RHIC and SPS energies. 

Now we summarize MICOR results on transverse slopes at projectile energies 
Ebeam = 40, 80, 158 AGcV at SPS. Figure^lshows the slopes in case of neglecting 
final state interactions. Figure [S] displays the slopes for including final state effects. 
In the latter case experimental data can be reproduced by transverse flow = 
0.38, 0.40, 0.44, respectively. These vt values are smaller than in the interactionless 
case, the difference is v^°^^^'^ —v'^^^^ ~ 0.04. It means that final state interactions 
increase transverse flow effect. Thus we can start from smaller quark flow values for 
the expanding quark matter, after switching on hadronic interactions the hadronic 
transverse momentum spectra will be reproduced. 

Inverse slopes on Figure |31 are fitted at low pT [pt < 1 GeV) . There is a 
small difference between data and calculations at smaller collision energies, even 
for the best choice of the vt values. Pion and kaon inverse slopes from MICOR 
are significantly higher than data at Ebeam = 40, 80 AGeV. At larger energy, 
Ebeam = 158 AGcV, the agreement is recovered using vt — 0.44 ± 0.02. 
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Fig. 6. The energy dependence of the average transverse flow velocities of the 
expanding quark matter in the MICOR model including final state effects (full 
boxes) and neglecting them (open boxes). 

Now the results of Figures [21 El and [3 can be combined into one figure to display 
energy dependence of transverse quark flow, vriy^SNN), see Figure O Including 
final state interactions (full boxes) we obtain smaller values for the quark flow than 
neglecting them in a fast hadronization scenario (open boxes). The uncertainties 
on the transverse flow values allow us to construct two interpretations. 

First we can assume a simple logarithmic energy dependence (dotted lines). If 
we neglect final state effects, then we can draw the following line to guide the eyes 

vt{V^) = 0.32 + 0.050 • ln(yi]^) . (4) 

Including final state effects the energy dependence of the quark transverse flow is 

vriV^) = 0.30 + 0.045 • ln(Vi]^) . (5) 

Here the energy y/sJjN is measured in GeV. 

On the other hand, our results obtained in discrete energy points may give a 
hint for a plateau between the largest CERN SPS and the smallest accessible RHIC 
energies (solid lines). Further studies are needed to decide if this hypothetical 
plateau exists at all and what are the theoretical consequences. The existence of 
such a plateau could yield important information about the equation of state of 
deconfined matter. 

The increasing transverse quark flow can not violate causality, thus it will sat- 
urate. However, this saturation seems to appear well beyond RHIC energies. 
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3. Conclusions 

In this paper we summarized our results on the energy dependence of quark trans- 
verse flow. Wc assumed that dcconfincd matter was produced at SPS and RHIC 
energies in central Pb+Pb and Au+Au collisions and an initial quark transverse 
flow has been developed before hadronization. Using the MICOR hadronization 
model wc could extract quark transverse flow parameters at phase transition. The 
comparison of calculated hadronic transverse momentum spectra and measured ones 
yields quark flow values with high precision in a wide energy region. Assuming fast 
hadronization and neglecting final state hadronic collisions, one can determine the 
wanted energy dependence for quark flow. We are pleased to see the good agree- 
ment between data and our quark based hadronization model calculations at low-pT 
{Pt < 1 GcV). 

Using the GROMIT transport code, one can investigate the influence of final 
state hadronic effects. In this case smaller quark flow values can be used to repro- 
duce experimental data. At first wc reproduce the proton transverse momentum 
spectra. Applying the extracted quark fiow, the pion and kaon spectra will be 
smaller than the experimental ones. However, this deviation allows us to include 
meson production from independent jet fragmentation, determined by pcrturbative 
QCD calculations. These perturbative contributions become comparable to the 
quark coalescence yield in the region > 1 GeV for pions and kaons. A detailed 
pcrturbative QCD calculation is needed to determine the interplay between pcrtur- 
bative, quark coalescence and hadronic contributions for light mesons. For heavy 
baryons and antibaryons the perturbative contributions are negligible in the region 
Pt < 4: GeV, where coalescence and secondary hadronic interactions dominate their 
momentum spectra. 

The energy dependence of the obtained quark transverse flow values can be 
described by a simple logarithmic increase. On the other hand, the high precision of 
the extracted flow values may indicate the possible existence of a plateau between 
SPS and available RHIC energies. Further analysis of the existing experimental 
data is needed to confirm this plateau. In parallel, new data would be needed in 
the energy range 20 < ^/snn < 130 GeV to determine the energy dependence more 
precisely. 
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